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Parity-time-symmetric ('PT-symmetric) optical waveguide couplers offer a great potential for 
future applications in integrated optics. Studies of nonlinear PT-symmetric couplers present new 
possibilities for ultracompact configurable all-optical signal processing. Here, we predict nonlinearly 
triggered transition from a full to a broken PT-symmetric regime in finite-size systems described 
by smooth permittivity profiles and, in particular, in a conventional discrete waveguide directional 
coupler configuration with a rectangular permittivity profile. These results suggest a practical route 
for experimental realization of such systems. 


One of the fundamental assumptions of quantum me¬ 
chanics states that operators corresponding to physical 
observables possess a real spectrum of eigenvalues [T]. 
The most common class of operators having purely real 
eigenvalues are Hermitian operators. However, in 1998 
Bender et al. showed that a more general class of 

non-Hermitian parity-time-symmetric (PT-symmetric) 
operators also possesses a real spectrum. Consequently, 
Hamiltonians with complex potentials V that fulfill the 
condition V{r) = ]/*(—r) have real eigenvalues and de¬ 
scribe physical phenomena. While the realization of such 
complex potentials in quantum mechanics is challenging, 
it is much easier to control them in optics. Indeed, the 
light propagation can be described by an optical analogue 
of the Schrodinger equation, where the complex dielectric 
permittivity distribution e(r) plays the role of potential 
and the imaginary part of the permittivity corresponds 
to gain or loss nig. 

The optical PT symmetry was experimentally demon¬ 
strated in linear optical waveguide couplers [6], which 
are important components for future fast, ultracompact, 
and configurable all-optical signal processing. Optical 
couplers with gain and loss were also studied in the non¬ 
linear regime [7] , where unidirectionality [8] and suppres¬ 
sion of time reversal [9] were shown. Phenomena related 
to the nonlinear PT symmetry were also studied in pe¬ 
riodic systems, including solitons [lOHISl, breathers m 
and their stability 

Optical PT-symmetric systems, where e(r) = e*(—r), 
undergo a transition with the change of the ratio be¬ 
tween the imaginary ejM and the real part eRE of the 
permittivity modulation depth (eiM/cRE)- Below a cer¬ 
tain threshold value of this ratio, the system is in the 
full PT-symmetric regime and has purely real eigenval¬ 
ues mnn]. For the eiu /^re ratio above this threshold 
value, the system is in the broken PT-symmetric regime. 
In this regime, a pair of modes possesses effective indices 
that are a pair of complex conjugate numbers, such that 
one mode experiences gain and the other loss. In lin¬ 
ear systems, the transition from the full to the broken 
PT-symmetric regime (hereafter simply PT transition) 
is usually controlled by changing the amount of gain and 



FIG. I. Geometry of the studied structures with their pa¬ 
rameters for different types of PT symmetric couplers. The 
size of rectangular waveguides is denoted by d, and the dis¬ 
tance between them is L. The size (period) of the cosine-like 
dimers is denoted by T. cb denotes the background relative 
permittivity; cre and cim denote the modulation amplitude 
of the real and the imaginary parts of relative permittivity, 
respectively. 

loss in the system (varying cim)- Nevertheless, the ratio 
eiM /cre is also influenced by the amplitude of cre , which 
in nonlinear systems can be controlled by varying the in¬ 
cident light intensity. Recently, it has been shown that 
the nonlinearity can trigger a PT transition in an infinite 
periodic array of PT-symmetric waveguides described by 
cosine-like permittivity distribution m- 

In this paper, we show that the nonlinear PT tran¬ 
sition can be observed in ultracompact nonlinear cou¬ 
plers consisting of one or two PT-symmetric waveguides, 
which can be readily integrated on a chip. Figureshows 
the three geometries under investigation: PT-symmetric 
couplers built of rectangular (a) and sinusoidal (b) waveg¬ 
uides; and a dimer (c), where both gain and loss are lo¬ 
cated in one waveguiding structure. 

To study the nonlinear light propagation in one¬ 
dimensional (ID) PT-symmetric structures with Kerr- 
type nonlinearity, we use the scalar wave equation for 
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the electric field E{x, z): 

+ k^eix) + k^aix)\Ef] E = 0, (1) 

where the operator jdx"^ +3“^ jdz^ denotes the 2D 

Laplacian operator, ko = 27r/A is the free-space wavevec- 
tor, and A is the free-space wavelength of light. The light 
propagates along the ^-direction, and both the structure 
and the field distributions are assumed to be invariant 
along the ^-direction. The linear relative complex per¬ 
mittivity distribution along the x-coordinate is described 
hy €{x) = €b -\- Ae(x), where denotes the background 
relative permittivity and Ae(x) describes the linear com¬ 
plex permittivity modulation depth. The Kerr nonlinear¬ 
ity strength is quantified by a{x) = eo5Re{e(x)}cn2(x), 
where eo denotes the vacuum permittivity, c denotes the 
speed of light, and n 2 is nonlinear parameter correspond¬ 
ing to cubic nonlinearity used in the definition of the in¬ 
tensity dependent refractive index Atinl = ^ 2 ! [21]. The 
light intensity is related to the electric field in the follow¬ 
ing way: I = eonBc\E\‘^/2^ where ub = yEs- The non¬ 
linearity a{x) in the structures studied here is nonzero 
only inside of the waveguides [i.e., if e{x) = e^, then 
a{x) = 0]. 

Using the slowly varying envelope approximation 
E{x,z) = ip{x, , Eq. Q is transformed into 

the (1 -h 1)D nonlinear Schrodinger equation 
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Equation is solved using the split-step Eourier method 
[22[ [23] in order to analyze the nonlinear dynamics of the 
light propagation in our T^T-symmetric structures. 

While some of the previous theoretical studies focused 
on the nonlinear light propagation in an infinite periodic 
cosine-like T^T-symmetric waveguide array [20], here we 
show that the rich variety of nonlinear phenomena can be 
obtained in simpler finite-size systems of waveguides that 
are more feasible from the laboratory viewpoint. Eirst, 
we study a nonlinear coupler built of two rectangular 
waveguides presented in Eig.j^a) and look for the nonlin¬ 
ear transition from the full to the broken T^T-symmetric 
regime. The structure with the following parameters is 
studied: d = 1 /im, L = 0.1 /im, = 2, and eRE = 0.02 
at the wavelength A = 0.63 /im. 

In order to study the properties of light in our VT- 
symmetric structures, we solve Eq. 0 using two different 
methods. Eirst, we study modal properties of the system 
in the linear regime [a{x) = 0]. We look for the field pro¬ 
files that propagate with an effective index ngff without 
changing their shape in the form E(x^ z) = 

This allows us to transform Eq. 0 into an eigenvalue 
problem for the modes of the waveguide described by 
e{x) and find the corresponding dispersion relations. We 
use these dispersion diagrams [shown in Eigs. ia), (b)] 
to locate the symmetry breaking threshold, where two 





FIG. 2. (a), (b) Linear dispersion relations for a pair of rectan¬ 
gular waveguides shown in Fig. [^a). The real (blue curves) 
and imaginary (red curves) parts of the effective index are 
shown as a function of the (a) imaginary eiM and (b) real 
cre part of the permittivity modulation depth. For (a) cre is 
fixed and equal to 0.02, whereas for (b) eiM is fixed end equal 
to 0.015. The black dashed line denotes ub- (c)-(f) Non¬ 
linear dynamics of light propagating in dimers built of (c)- 
(e) two rectangular waveguides (where 77-2 = 10“^® m^/W) 
and (f) two cosine-like waveguides (where a — 10“^^ m^/V). 
(c), (e), (f) The intensity distribution I[x^z) and (d) the evo¬ 
lution of the x-coordinate of the center of mass. The geome¬ 
try of the waveguides is presented in Fig.[^a), (b) and in the 
figure insets. The initial power density Pq is set to (c), (d) 
10^° W/m, (e) 3 • 10^° W/m, and (f) 5 • 10® W/m. White 
horizontal lines indicate waveguide boundaries of rectangular 
waveguides. 


modes with purely real effective indices are transformed 
into a pair of modes with complex conjugate effective 
indices. 

For the coupler parameters given above, this VT tran¬ 
sition occurs at cim = 0.0149. Here, we study the sys¬ 
tem slightly above the phase transition (at cjm = 0.015, 
which corresponds to the broken T^T-symmetric regime 
for low light intensity). The nonlinear dynamics shown 
in Figs.j^c), (d) proves that it is possible to observe the 
nonlinearly triggered VT transition in the dimer built 
of two rectangular waveguides [(shown in Fig. [^a)]. In 
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this and all of the other configurations studied here, 
we excite the gain mode that initially increases its en¬ 
ergy. At the beginning of the propagation, the light in¬ 
tensity is low and the system is located in the broken 
T^T-symmetric regime. The light is mostly located in 
the gain waveguide, as it can be seen in Fig. [^d) show¬ 
ing the x-coordinate of the center of mass (COM) calcu¬ 
lated as xl{x)dx/ /^^/(x)dx. With the increase 
of the light intensity, the real part of the permittivity 
of the waveguides also increases. Since the VT transi¬ 
tion depends on the ratio between the real and imaginary 
parts of permittivities [20] , the increase of eRE drives the 
system back to the full T^T-symmetric regime [the dis¬ 
persion relation neff(eRE) presented in Fig. [^b) shows 
that with the increase of eRE, the system is transformed 
from the broken to the full VT -symmetric regime]. After 
the VT transition occurs, the light intensity in the gain 
waveguide reaches its maximum, and then it starts to di¬ 
minish. At the same time, energy is coupled to the lossy 
waveguide in which the maximum intensity is reached 
slightly later and then also starts to decay. With the de¬ 
crease of the total intensity in the coupler, accompanied 
by the decrease of the nonlinear permittivity modula¬ 
tion depth, the system is transformed back to the broken 
T^T-symmetric regime. This cycle repeats itself during 
the propagation. 

Figure [^e) presents the light evolution in the same 
system as the one in Figs, [^a), (b) but for a higher ex¬ 
citation power density Pq = I{x^z = 0)dx. We can 
observe that with the increase of the initial power, the 
period of oscillations decreases and the maximum light 
intensity increases. At the low power level, the oscillation 
period is approximately equal to 1.5 mm, and the peak 
power corresponds to the nonlinear permittivity modula¬ 
tion depth of the order of 1.7 • 10“^. For high power, the 
period decreases to 0.6 mm, and the maximum nonlinear 
permittivity modulation depth increases to 3.5 • 10“^. 

Let us now compare these results for the discrete rect¬ 
angular waveguide-based coupler with those for the cou¬ 
pler with a smooth cosine-like profile described by 
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FIG. 3. Nonlinear dynamics of light propagating in a single 
dimer described by Eq. (a), (c) The intensity distribution 
/(x,z), (b) the evolution of the x-coordinate of the center 
of mass corresponding to subplot (a). The geometry of the 
structure is shown in Fig. [^c) and in the inset of supblot 
(a). The parameters are: cre = 0.05, a — 10“^^ m^/V^, 
initial power density Po = 10® W/m, and the imaginary part 
of permittivity eiM equals (a), (b) 0.031 and (c) 0.035. 


the nonlinearity-induced VT transition is not sensitive to 
the exact refractive index profile of the structure. These 
results are likely to be important for the experimental 
observation of the predicted phenomena. 

Now, having established the similarity between light 
propagation in PT-symmetric structures defined by rect¬ 
angular and smooth permittivity profiles, we consider the 
effect of different distributions of both the real and imag¬ 
inary parts of the dielectric permittivity on the VT tran¬ 
sition dynamics. In particular, we investigate a finite-size 
PT-symmetric dimer that can be described by 


Ae(x) = cre cos^(2flx) + iejM sin[2fl(x — 7r/4)] (3) 


Ae(x) = cre cos^(flx) + ieiM sin(2flx), (4) 


and shown in Fig. [^b). Here, = tt/F, F denotes the 
full width of the dimer, and both cre and cim are real 
quantities. While such smooth cosine-like profiles are 
often used in theoretical studies of periodic structures 
and waveguides, it maybe challenging to realize such pro¬ 
files in laboratory experiments. The cosine-like coupler 
parameters are: cb = 2, cre = 0.03, and F = 3 /im. 
The imaginary part of the permittivity is chosen to be 
Gm = 0.0162, which is just above the VT transition point 
located at ejM = 0.016. It is noteworthy, that the dy¬ 
namics of nonlinear wave propagation in such a struc¬ 
ture, shown in Fig. |^f), resembles that of the discrete 
waveguides based coupler [see Fig. [^e)], suggesting that 


with Kerr-type nonlinearity. The geometry of a single 
dimer described by Eq. @ is shown in Fig.j^c), and the 
parameters are = 2, eRE = 0.05, and F = 3 /rm. Here, 
the transition from the full to the broken PT-symmetric 
regime is also observed. The threshold in the linear case 
is located at ejM = 0.0305. 

The nonlinear propagation of light in a single dimer 
is presented in Figs. ia), (b). It resembles the prop¬ 
agation of light in the unit cell of a periodic array of 
dimers presented in Figs. 3(d), (f) in Ref. [20]. However, 
in Ref. [20] the parameters of the array were chosen in 
such a way that the energy transfer from the gain region 
to the lossy one was visible only from the plot showing 
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the center of mass and not from the color map present¬ 
ing the intensity evolution during the propagation. On 
the contrary, for the parameters chosen here, this energy 
transfer, illustrated as the displacement of the center of 
mass in Fig.|^b), is also clearly visible from the intensity 
distribution presented in Fig. [^a). Indeed, the intensity 
distribution pattern that changes periodically along the 
2 ;-coordinate is visibly tilted away from the vertical axis, 
indicating the energy transfer from the top part of the 
dimer (gain region) to the bottom part (loss region). 

Figures [^a), (b) show that the system is initially (for 
low light intensity) in the broken T^T-symmetric regime. 
The gain mode of the system is excited, and the energy 
of this mode increases. Similar to the case shown in 
Fig. ID this increase is accompanied by the increase of 
the real part of the permittivity due to the Kerr effect. 
When the permittivity modulation depth is sufficiently 
high, it brings the system back to the full T^T-symmetric 
regime. Initially, the energy is stored in the gain mode 
of the broken T^T-symmetric system. In the full VT- 
symmetric regime, this mode is no longer an eigenmode 
of the system, and the distribution of light (still experi¬ 
encing gain as it is localized mostly in the gain region) 
is attracted to the center of the dimer, where the per¬ 
mittivity is the highest [20]. Therefore, the beam pos¬ 
sesses a nonzero transverse momentum directed towards 
negative x-direction. Due to this momentum, the beam 
crosses the center of the dimer and as a result is mostly 
localized in the region with loss. Consequently, the to¬ 
tal energy rapidly decreases and the system is trans¬ 
formed back to the broken T^T-symmetric regime. In 
this regime, the part of the energy stored in the loss 
mode vanishes quickly and the part stored in the gain 
mode starts to grow again. This cycle repeats, as it is 
seen in Figs. a), (b). 

It is noteworthy that in the case of a single dimer, for 
which the parameters are chosen in such a way, that the 
system is higher above the VT transition (ejM = 0.035 
while it was 0.031 in the example discussed above), the 
nonlinear transition also occurs, but it has a different 
character. The propagation of light in such a system is 
presented in Fig. [^c). We see that in the initial phase 
of the propagation (up to 0.3 mm), the oscillations are 
aperiodic, and the period of oscillations decreases. With 
the decrease of the oscillation period, the energy is trans¬ 
ferred from the gain part of the dimer to the lossy part 
within a shorter distance. At the beginning of the prop¬ 
agation, the transverse component (toward negative x 
values) of the beam momentum is low enough so that 
the beam stays confined in the dimer. During the propa¬ 
gation, this momentum increases and at the propagation 
distance x ^ 0.3 mm, it causes a significant quantity of 
energy to be outcoupled from the dimer and to propa¬ 
gate toward negative values of x. After this point, the 
propagation pattern becomes more regular: the period of 
oscillations is kept constant, and only the peak intensity 


reached during the oscillations changes (the variations of 
the peak intensity are at the level of 50%). The irregular¬ 
ity of the propagation in the system located higher above 
the VT transition can be explained by the fact that the 
light intensity here is 20 times larger than that for the 
system closer to the VT transition. At the same time, 
the period of the oscillations is about 4 times smaller, and 
consequently, the changes in the light intensity are much 
more abrupt than in the system closer to the VT transi¬ 
tion. These rapid changes cause more irregular behavior 
of the system high above the VT transition. 

Here, we briefly compare the field distribution in a pair 
of waveguides shown in Fig. [^and the single cosine-like 
dimer shown in Fig.|^ The main difference lies in the fact 
that in the case of the single cosine-like dimer, the field 
possesses a single intensity maximum whose location is 
described by the center of mass that shifts from the gain 
to the loss part of the waveguide, whereas in the pair of 
waveguides the field distribution has two maxima—one 
in each waveguide. In the system of two waveguides, it 
is the ratio between these maximum intensities (or more 
precisely the integrated power in each of the waveguides) 
that determines whether gain or loss dominates at a given 
propagation distance. 

In conclusions, we have studied nonlinear dynamics 
of light propagating in finite-size parity-time-symmetric 
optical couplers. We have shown that the nonlinearly 
triggered transition from the full to the broken VT- 
symmetric regime can be observed for various finite-size 
couplers built of two rectangular and cosine-like waveg¬ 
uides. We have shown that the nonlinear ly triggered 
transition between the full and the broken T^T-symmetric 
regime also occurs in a single T^T-symmetric cosine-like 
dimer. Moreover, we have found that the light propaga¬ 
tion in smooth and discrete waveguide couplers is qualita¬ 
tively similar. These results suggest a practical and sim¬ 
ple route to experimental verification and applications of 
this kind of nonlinear VT transition. 
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